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ABSTRACT 


In the present work, a methodology for Elastodynamic 
analysis of the valve train mechanism has been presented. 

The valve train mechanism is essentially ,a cam-nodulated 
linkage. This mechonism is widely used in I.C, engines and 
forms an important subsystem, A model for the elastodyna- 
mic analysis of a higher pair has been proposed . Using 
this proposed model for tho higher pair and the conventional 
models for lower pairs of the valve train mechanism, it has 
been shown as to how to obtain the error between the cam 
command and follower response. The error is a function of 
the input speed and this dependence has been studied, A 
computer program based on the methodology proposed has been 
developed, tested. Using this programme, it is possible to 
study the influence of various design parameter on the 
performance of the valve train mechanism. 
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CHAPTER I 
INTRODUCTION 

1.1 ELASTODYNAf^ilCS OF PLANOil MECHANISM} 

The kinematic analysis and synthesis of 
any planac mechanism has become a powerful design tool for 
engineers. The technique of kinematic analysis, hov/ever, 
suffers from a major drawback* This is the ’rigidity’ 
assumption which prevails throughout the literature with 
few notable exceptions. Mechanisms consisting of gears, 
links, sliders etc, in real world application are not at 
all rigid* They are sensitive to clastic deflections when 
they are subjected to higl static or dynamic forces. The 
dynamic forces play an insignificant role when the machine 
or mechanism is operated at low speed range* Consequently, 
a designer does not need to concern himself about the inherent 
elasticity present within the system. The idealization of 
the mechanism as consisting of ’purely rigid members’ in the 
analysis or design process does yield results v/ithin the 
designers satisfactory level as long as the operating speed 
IS kept low and so long the accuracy required in the rela- 
tionship of input and output is not very high. 
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The situation turns out to be significantly 
different in the case when the same mechanism is operated 
at higher speeds. Dynamic forces increase with the increase 
of acceleration of different members. For sufficiently larger 
speeds, these force may even play a predominant role. Thus 
the rigidity assumption, generates inaccurate analysis and 
synthesis which may lead to improper functioning of the 
machine and eventually failure of the members particularly 
at high speeds. A machine or mechanism which is kinematic- 
ally adequate at slow speeds may have sufficiently differemc 
motion at high speed so that it can no longer perform its 
intended function. Apart from this, in high speed mechanisms, 
the chances of the mechanism being subjected to resonant 
condition is also quite high* 

V/ithin the domain of present day automation, 
machineries are required to o’^erate at higher speeds and at 
axnimum power. The accuracy requirements in such machines are 
also showing an increasing trend. Elastic deflections which 
are the consequences of increased speeds and loads may cause 
inaccuracies of position in addition to noise and fatigue. 

From the above facts, it becomes thus apparent 
that the usual design or analysis procedures may not yield 
satisfactory results. Need for newer methods and advanced 
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techniques was felb for analysis aid design of mechanisms 
Vi/ith elastic members. Literature has already recognized a 
considerable amount of research work in this area. The 
outcome of these efforts has maoo a way to a newer and 
advanced method, Kineto~elas todynamics , an analysis and 
synthesis technique which can take care of kinematics as well 
as finite rigidity of the members. It has become an efficient 
and powerful bool to the designers, 

Kineto-elas todynamics in general, is the exami- 
nation of displacement, velocities, accelerations, stresses, 
strains etc* of a moving elastic mechanism, effects of 
elastic deformations upon the inertia forces are included in 
the analvsis* Whereas, the Elastoc'vnamic analysis of a 
mechanism defines the same thing as the previous one except 
that the inertia forces are calculated by assuming the 
members rigid. Both the methods take into account the mass 
distribution of the links. The principles and concepts of 
finite element method used in structural dynamic analysis 
IS greatly utilised in performing the kineto-elastodynamic 
analysis of any planar mechanism. 

Over the past decade, since its origin, lot of 
improvements have been carried out towards the generalisation 
of KED analysis technique for studying lower pair planar 



mechanisms 


Most of the researchers in this field have 


developed methods or improved over the previous methods 
specific application to the planer four bar mechanisms 
(e.g., 4-bar slider-crank mechanism, crank-rocker mechanism 
4-bar path generating mechanism etc.). Their applications 
were limited to planer lower-pair mechanism. On the other 
hand, higher-pair mechanisms did not get much attention 
on the employment of the analysis method. 

Cam-modulated mechanisms are noteworthy examples 
of higher-pair mechanism. In the construction of modern 
machines their importance are obvious. They find their use 
immensely in automdtive engines, aircraft engines, textile 
machineries, automatic machines etc. Overhead valve-train 
linkage used in automative engine is one of tho simplest 
example of a cam-modulated linkage. In the context of smooth 
operation, efficient performance of the engine tho valve- 
train mechanism^ gets due importance. The overall performance 
of mechanism depends largly on tho proper design of the cam 
profile. Hence, a comprehensive and accurate analysis comp- 
ared to the previously used ones of the complete mechanism 
is essential which helps in acheiving better performance. 
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1.2 OVERHEAD VALVE- TRAIN IN AUTOMOTIVE ENGINE : 

A typical cam-driven overhead valve-train linkage 
used in automative engine is shown in Fig. 1.1. This mecha- 
nism consists of a pushrod, rocker arm, valve, value-spring 
apart from the cam and follower. The movement of the valve 
IS governed by the cam profile. The valve controls the flow 
of air fuel mixture into the engine cylinder. The amount of 
air fuel mixture getting into the cvlinder determines to some 
extent the thermodynamic efficiencv of the engine, thus the 
importance of the proper valve movement. 

In almost all the automotive engines the response 
of the follovwer end that is of the valve end does not remain 
same as that of machined on the cam profile, ( i.e., cam 
command). This is inevitable particularly at higher speed. 

The discrepancy in cam command and follower response is a 
result of the high flexibility or elasticity of the follower 
train and the magnitude of dvnamic forces prevalent in the 
mechanism as a consenuenco of higher speed. In the context 
of accurate design of cam profile correct estimation of the 
difference between cam command and resnonse is essential. 

Thus a need is felt to in^^estigate the discre- 
pancy through comprehensive analvsis utilising the newer 
methods, like elastodvnamics , 
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1,3 LITERATURE REVIE'Us 

During tho earlier stage the works that were 
conductecU the researcher analysec’ the mechanisms based on 
the mathematical model with only one elastic member [lO]. 

For instance, in 1965, Neubauer, Cohen and Hall [lO] analy- 
sed a slider-crank mechanism containing an elastic coupler 
rod. In the next decade the analvsis method developed, 
assumed either most or all the momoers were elastic. 

Vi/infrey [21] performed an analysis of longitu- 
dinal, transverse and torsional displacements of a general 
mechanism. He v^as the first to use finite element method 
in the area of mechanism analysis. For each planer link he 
defined 6 clastic cordinates ( 3 at each end of the link ) 
and 3 rigid body coordinates. 

Erdman, Sandor and Oakberg [s] developed a gener' 
al method for kineto-elastodynamic analvsis and synthesis 
of planor mechanism. They incorporated flexibility approach 
of the finite element method in the analysis. The input 
crank is considered to behave like a cantilever beam thus 
converting the mechanism to ^ structure ’ and a system 
flexibility matrix is developed, A Kineto-elastodynamic 
stretch rotation operator ( KEDSRO ) describing the elastic 
motion of a single point is used. This finds out the 
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deflections by an iterative procedure, which starts from 
rigid body inertia forces. 

In 1973, Iman,Sandor and Kramer [7] came up with 
a method of deflection analysis of general planoir single and 
multi— loop linkages. The extension of permutation vector 
IS used to carry out the assembly of elements. This method 
numbers the system co-ordinates automatically and it d-oes 
not offer any special advantages in the sense that system 
coordinates are less in number in case of a mechanism. They 
also use for the fjrst time the rate of change of eigenvalue 
and eigenvectors in mechanism analysis in an effort to reduce 
the computational time, A reduction by a factor of 3 is 
found by them. The force vector is calculated bv multiplying 
the system mass matrix with the negative of the acceleration 
( rigid body ) vector and assumed to remain constant during 
each interval. Input crank is considered as a cantilever 
beam since no system co-ordinate is defined at the support 
end. A method of computing dynamic stresses is also discu- 
ssed in this paper, 

Mdha , Erdman and Frohrib [8] described general 
procedures for studying high speed elastic linkages using 
finite element displacement method of structural analysis 
technique. The expressions for nodal displacements. 
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accelerations etc. include the terns coupling the rigid~body 
and elastic motions. The effects of additional acceleration 
terms are neglected in the final form of these expressions. 

The input crank here also is modeled as a cantilever beam, 

Winfrey [l9j reported a method for dynamic 
analysis of flexible mechanism by reduction of co-ordinates. 
The total number of svstera co-ordinates is reduced to only 
one cordinate at the point of design interest. Thus a consi- 
derable amount of computational time is saved by solving 
a single degree of freedom system of the same elastic mecha- 
nism, Results he obtained'’, as expected differs from that 
of the complete solution. 

In his analysis method Nath [9] treated each 
link with no, of subdivisions for accurate modelling as 
compared to the previous works. A new and efficient technique 
is developed to eliminate the rigid body degrees of freedom 
instead of taking the input link as a cantilever beam. V\fhile 
deriving the element property matrices he takes into account 
the effect of rotary inertia and shear. To construct the 
svstem matrices code-system of assembly is adopted to incor- 
porate various mechanism conditions. 

Erdman and Sandor [6] presented a review of the 
state of the art and trends in kineto-elastodynamics. The 
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nomenclature used in KED is consolidated in this v;ork in 
organised form. 

Winfrey, Anderson and Gnilka §-c>j studied an 
elastic machinery with clearance. The analysis approach 
is presented with elastic model citing the example of a 
cam-driven autometive overhead valve- train linkage. The 
elastic model of the follower train is formulated using the 
finite element technique for structural analysis. The entire 
analysis process assumes constant geometry of the mechanism. 

He models the component parts as simple prismatic elements. 

The proDlem of impact is also investigated in tnis work taking 
help of the elastic model prepared for the follwer train. 

Previous works as discussed in the above-mentio- 
ned literature contributed towards the establishment of the 
foundation of the kineto-elastodynamic technique. Their 
applications are confined around primitive four-bar mechanism 
and also to multi-loop linkages in few instances. From the 
point of view of applying the kineto-elastodynamic analysis 
to cam-activated mechanism almost no prior effort is made. 

Of course, works have been carried out in the past decades 
by numerous researchers to model a cam-follower system. 

Many works have been published on the study of dynamic 
response of different types cam-follower system with the 
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aid of conventional single-degree or multiple degree of 
freedom mass, spring, damper representation. 

Barken [2] dynamically modeled an automative 
overhead valve train using the lumped mass approach. Ini- 
tially he forms a multidegree of freedom model containing 
massless springs and rigid masses. Then he reduces it to a 
single-degree of freedom model and solves to find dvnamic 
response of the svstem . Stoddart [l4] introduced a new 
method, namely ( Polydyne cam design ) to design a cam pro- 
file for an automotive engine accounting the elasticity of 
the follower train. In the present work an attempt 3S made 
to estimate the elastic deflection of a cam-actuated linkage 
exploiting elastodynamic technique, 

1,4 OBJECTIVE OF THE PRESE-.T h'ORK ; 

The analysis or synthesis of a cam-follower 
system performed by numerous workers in the previous years 
utilizing lumped mass approach and mass-spring damper model 
suffers from few drawbacks. These can bo stated as lack of 
generality in approach, inaccuracy incurred in calculating 
deflections etc* The concept of finite element method with 
its versatility may be used to eliminate the above shortco- 
mings in the context of analysis of cam mechanism. 
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Analyzing a cam-actuated follov/er train with 
the aid of elastodynamic technique possess certain difficu- 
lties. Being a higher pair the cam-follower cannot directly 
be ficted into the me thodologv of elastodynamic modelling 
technique for the purpose of modelling them# 

Unlike the links in a lower-pair mechanism the 
kno\m structural beam elements cannot possibly be applied 
to model the cam-follovver pair directly* Obviously* it 
becomes necessary to obtain an appropriate equivalent model 
of the cam- follower pair which When combined With the rest 
of tne follower linkage enables one to exploit the potenti- 
ality of elasto-dynamic analysis technique and eventually 
to perform the analysis. 

The objective of the present work may be narrated 

as follows; 

(i) an attempt will be made to represent a cam- 
•Pollo'fer pair to its corresponding lower-pair elements which 
can approximately ( if not accurately ) becomes equivalent as 
far as the elastic and dynamic behaviour of the actual system 
is concerned. The lower-pair equivalent model will also be 
made kinematically equivalent. It becomes apparent that the 
finite rigidity of the cam body comes into picture and in a 
way forms a part of total system stiffness. 
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(ii) to formulate tho system model interms of system 
co-ordinates in order to perform the deflection anal''''sis, Tne 
elastic behaviour of the system ’will oe examined at various 
operating speed. 

The analyses will be carried out for an automoti- 
ve overhead valve-train linkage which is a cam-modulated higher 
pair mechanism, A radial cam with translating roller I'ollower 
is chooser for the sake of simplicity in illustrating the 
analvsis procedure throughout. 

1.5 LIIdJTATICNS: 

The analyses aoproach tnat is follovved has got 
few limitations and assumptions. 

In formulatjng the system model For the complete 
valve train the input crank in the equivalent mechanism is 
modeled as a cantilever beam. This bypasses the purpose of 
eliminating the rigid body degrees of freedom and the 
mechanism is not analysed under actual boundary conditions, 

Tho force vector in solving the system equations 
of motion consists of the contribution from inertia forces 
arising due to gross or rigid body motion of the mechanism in 
the course of motion. This does not take into account the 
actual distribution of the rigid body inertia forces over the 
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lengths of the element. Any external forces if present, is 
also included in Torce vector. 

The assumptions that are taken are : 

(i) the input angular velocity is assumed to remain 

constant. 

(ii) the rigid body motion and elastic vibrations are 
seperable , 

(iii) the deflections are small so that the linear 
theory can be applied throughout. 

(iv)the effects of clearances, tolerances are 


neglected 
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CHAPTER II 

ELASTODYNAMIC MODELLING OF VALVE TRAIN 

2.1 KINEMATICS : 

As indicated in the Chapter I the type of 
cam-follower pair that is being considered in the present 
work a radial cam with translating roller follower. In order 
to carry out the elasto-dynamic analysis of the entire 
mechanism as a first step it is necessary to model the cam- 
follower pair first. Most of the lower pair mechanisms 
(like, slider-crank, crank-rocker etc.) are comprised of 
binary links mainly and of ternary links in some cases. 
Analysing those mechanisms, a binary link can very simply 
be represented by a structural beam element. The ternary 
link can also be represented by attaching appropriate degrees 
of freedom at the nodes. But a cam-follower pair which is 
a higher-pair ( that means the nature of contact between 
them is either point or st.line ) cannot directly be modeled 
using simple beam elements. Cam is basically a disc and follo- 
wers also do appear in various forms. 

Thus difficulties are encountered in modell- 
ing the cam follower pair. But a higher pair can also be 
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converted to a lower pair such that the latter becomes 
instantaneously equivalent kinematically to the former one. 
That is at the point of contact the kinematic quantities 
like displacements, velocity, and accelerations are all 
identical instantaneously for both the higher-pair and 
lower pair. Similarly, a cam-roller follower can always 
be transformed into an equivalent lower pair mechanism to 
obtain a instantaneous kinematic equivalency at the point 
of contact. Having made the kinematic equivalency by 
converting to a lower-pair structure it can be possible 
now to start with analysis if the dynamic equivalency of 
the two svstem is achieved. In the following section of 
this chapter ( section 2,2 ) an attempt is made to obtain 
an equivalent elasto-dynamic model of the cam-pair at least 
approximately. 


Let us consider a radial cam and a transla- 
ting roller follower ( see Fig, 2.1 ) . The direction 
of the follower motion is conciding with the T'-axis of the 
cam, X-Y is the fixed axis and Xj^- y^^ is a moving co-ordinate 
frame embedded on the cam profile. The follower displacement 
is dictated by the characteristic shape of the cam profile . 
Let P be the point of contact at the instant shown when the 
Xj-yi frame inclines say, at some angle cj) with X-Y . To 
create an instantaneously equivalent lower-pair mechanism the 




A RADIAL CAM AND TRANSLATING 
ROLLER FOLLOWER . 


INSTANTANEOUS KINfMATIC 
EQUIVALENT LINKAGE 


FIG. 2.1 


FIG. 2.2 
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centre of curvature C of the cam~prof ile is located. For 
the roller follower the same is the point Q. 

It can be shown that [ll] OCPQ which represents 
a slider-crank mechanism ( see Fig, 2,lj2,2)is equivalent 
instantaneous 4 bar mechanism to the situation under consi- 
deration ,Thefollower motion and the slider motion are 
identical for the instant considered. Distance PQ for a roll- 
er follower always remains constant. The location C,thus the 
distance OC and CP varies with the cam profile, A harmonic 
cam with constant radius of curvature can be permanently 
represented by a slider-crank mechanism with constant link 
lengths, V/ith arbitrary cam profile tho link lengths OC and 
CQ changes when the actual cam rotate through any definite 
angle. 


2.1.1 RIGID BODY DISPLACEMENT : 

When a mechanism undergoes rigid body displacement 
no elasiticy or flexibility comes into picture. Members are 
assumed purely rigid. The functioning of the overhead valve 
train ( see Fig, 1,1) is as follows. The pushrod is constrained 
to move only along the guide. The rotation of the cam causes 
the pushrod move up and down depending on nature of the cam 
profile. During the rise period ( when the valve opens ) IL- 
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pushrod is raised up by the cam action. This causes one 
end of the rocker arm to move up. The valve side end of the 
rocker arm pushes down on the valve stem-causing the valve 
to move down or open. The valve takes the close position by 
the spring action when the cam passes through return period. 
The displacement at the cam end is faithfully transmitted to 
the volve end multiplied by a-f actor (lever ratio ), 

Let the displacement at the cam end be dc, 

measured from a datum ( say the base circle), when the cam 

has turned through an angle of ^ degree. At this instant dc 

is found from the profile curve, d can also bo obtained 

c 

from the equivalent linkage. The valve-end displacement is 
given as 


d^^ = d^ X c, (2.1) 

where c- rocker arm lever ratio 
d^- valve end displacement 
assuming small rigid body displacement. 

From the Fig, 2.2 , the following expression 
for d is written, valid for the instant shown 

dc = ^2 Sin T+ Sin 0 - 


( 2 . 2 ) 
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where, 

r 2 “ connecting rod length 
rj^ - crank length 

shown in Fig. 2.2 
G 

- base circle radius 

2.2 ELASTO-DYNAiMIC MODELLING OF THE CAiVi-FOLLO' 'ER PAIR ; 

In the previous section of this chapter a 
kinematically eguivalent lower-pair raechanism has been found 
out w/hich IS a slider crank 4-bar linkage, ]n general the two 
link lengths (Fig, 2,2 ) r^^ ( of the crank ) and r 2 ( of the 

connecting rod ) are expected to vary because the centre of 
curvature, is not going to remain fixed. Kinematic equiva- 
lency assumes links to be rigid. The aim of this section is 
to form a equivalent model of the cam-follower system for 
the purpose of elas to -dynamic analvsis. Basically, the task 
of modelling finally comes down to modelling a lower-pair 
system incorporating appropriately the elastic and dynamic 
behaviour of the actual svstem. 

It is observed (Fig. 2,1 ) that in the actual 
system the cam is purely rotating body, rotating about the 
axis passing through 0 and perpendicular to the plane of the 
paper. The crank in the equivalent slider-crank mechanism, is 
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purely rotating member. The ‘crank’ and the cam are revolving 
about the same centre of rotation. From the fact it can be 
assumed that the ’crank’ and the cam are maintaining rotational 
equivalency if their mass moment of inertia are equalised about 
the same centre of rotation. 

In other words, it is leant that the ’crank* is 
behaving like a ’rod like cam’. Equating mass moment of iner- 
tia of the two purely rotating body both are moving at same 
angular velocity, implies their dynamic equivalency. Thus the 
crank of the kinematically equivalent mechanism can have 
some mass. 


Let the mass moment of inertia of the actual cam 
about the centre of rotation 0 - and 

= mass moment of inertia of the crank about 0. 
cr 

Consider the crank is having uniform circular cross 
section area of radius r. If the cam profile is known comple- 
tely the it’s mass moment of inertia can be computed. For a 
thin disc/ plate of uniform thickness ’t’ and made up of a 
homogeneous material of density ’y* the polar moment of iner- 
tia I^^ ’ is given by [s], 

T = I ’ = V t J 

c cc ,mass ' * c,area 


.. .. ( 2 . 3 ) 
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where, 

'^c,area ~ polar moment of inertia of area of 
the plate about point 0. 


^c,area ® given cam profile can be computed 
numerically ( discussed in Appendix I ). Crank is considered 
as a circular cylinder. Mass moment of inertia of it about 
0 is expressed as [3 ] 


Icr = l2 ^ + ^^1 ) ^2.4) 

where , 

m = mass of the crank, 
r = radius of cylinder 

Equating ( 2,3) and ( 2,4 ) 

Ic = if ( 

= { 3 I + 4x2) 

Where , 

K = Tryrj^ 

2 

or, 12KI^= + AKxf m .....(2.5a) 

The equation (2, 5a ) is solved for unknown m as 

follows, 

-4r^ K + 4^(r^ K)^+9KI^ ,...(2.6) 


6 
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Feasibile value of 


-4r^K + 


9 


(rjK) 


+ 9KI, 


... (2,7) 


6 


Though in the actual system the body of cam has 
got very high stiffness compared to the other member in the 
mechanism. An estimation of this can be made approximately 
so that the 'crank' can be made to possessf inite amount of 
bending stiffness. The cam is assumed to be a cantilever of 
varying cross-section ( see Fig. 2. 3, 2. 4). Under the load P 
the maximum deflection cam be found as 


d = 


2 / 1 


[ ln(|) -(1- a ) + i ( 1 - )] 




( 2 . 8 ) 


(see appendix II) 


where , 


t = uniform thickness of the cam 


(_b-) 


From this expression the stiffness Kj^ is calculated 

as = I E.t/ -^^3[ In ( •^-) -(l-g) + | (1- ^ ) ] 

c b 

( 2.9) 

The connecting rod of the equivalent linkage is 
assumed massless. Since the actual cam and the crank are 
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made rotationally equivalent and the cam’s counterpart 
mass in the crank is assigned, the connecting roo can be 
taken as massless* But it is permitted to deflect elastica" 
lly only in its longitudinal direction. The stiffness of 
this link can be ettributed due to Hertzian contactness at 
the cam-follower interface. 


In actual cam-follower svstem a cylinder to 
cylinder contact 3 .s considered ( 'see irig. 2.5) [4]. 


The maximum normal stress 0, 


max 


is given as 


<^max= 


(Tl+Yg) 

Y1Y2 


1 - 


^1 


■u. 




1/2 


. . (2.10) 


-2 

where, = Young's modulus and Poissons ratio 

for the body 1 ( roller ) respectively 

Similarly £ 2 ?!-’'^ for the boc^y 2 { cam). 

Taking similar materials and |.i^=p 2 ='p 

expression (2.10) takes the form as in ( 2,11) 

1/2 


0 =,564 

max 


pt 

Y1+Y2' 

"liizr’ 


...( 2.11) 



26 


The stiffness Kj_j assigned to the connecting rod 
computed as 


K, 


H 


* • • « ( 2 » 12 ) 


.564 [- 


2(l-p^) 


,]l/2 


1 

fE 


Y = 


Vl2 

Yi Y2 


It is to be noted that for every geometric 
configuration of the cam the Ygj consequently Kj^ 
values are going to change. 

Thus, an estimation of the mass and stiffness 
of the two links in the original kinematically equivalent 
linkage could be made, and an approximate model for the 
actual cam follower system is constructed. This model will 
be utilised in formulating the system equations. 


2.3 ELASTODYNAfviIC MODELLING OF THE LINKAGE : 


Elastodynamic analysis is the examination of 
displacement, velocity, stresses etc, of a moving elastic 
mechanism and inertia forces are calculated by assuming all 
the members rigid. Finite element method is utilised to 
perform the analysis. Modelling of the entire mechanism 



27 


that is conceiving the mechanism as collection of elements 
is essential as a first step. It is mentioned here (again) that 
the input crank is modeled as cantilever beam which elimina- 
tes the need for removing rigid body degrees of freedom. Elas- 
todynamic analysis of a mechanism in motion helps in gaining 
insight into the elastic behaviour of different members 
under the dynamic conditions, 

2»3*1 ELEMENT MATRICES: 

In order to construct the finite element model 

each member of the linkage is taken as an element. Usually 

a link is modeled as a planer beam element used in structure, 

beam 

For a generalized planar £e lament 6 degrees of freedom, 3 at 

each end (node) is sufficient to completely describe elastic 
deformation of the element. This is shown in(Fig,2,6) in 
a 2-D frame, u^^ and u^ are axial elastic deformation coor- 
dinates. U 2 »U 2 and u^,u^ account for the transverse deflec- 
tion characteristic of the element Uj^ to u^ forms a genera- 
lised set of coordinates. They are defined as function of 
co-ordinates and time in the case of deriving element 
matrices for a link. The co-ordinates defining the axial 
deflections, a linear interpolation model ( mode shape ) is 
used. That is linear function of the co-ordinate defined 
along the element length. For the transverse displacements 
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a culDic interpolation function is assumed. Using linear 
superposition theory the tvjo deflection patterns are superi™ 
mposed so that complete elastic behaviour of the element is 
described. From equilibrium equations the element mass and 
stiffness metrics are derived [QJ* They are same with the 
stiffness and mass matrices of a planar frame element with 
6 d.o.f. 


The generalised form of those matrices are given 


below. 

It IS 

noted that 

these are 

element oriented element 

matrices and ; 

independent 

of svstem 

configuration. 


The e lement 

stiffness 1 

matrix [k] is 


EA/L 



— 


0 

12EI/L^ 




0 

6EI/L^ 

4EI/L 

Symmetric 

[k]= 

-ea/l 

0 

0 

EA/L 


0 

-12EI/L^ 

-6EI/L^ 

0 12EI/L^ 


0 

6EI/L^ 

2EI/L 

0 -6EI/L^ AEl/L 


(2.13) 
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The element mass matrix [m] is 


1/3 






0 

13/35 





0 

llL/210 

lVi0 5 




1/6 

0 

0 

1/3 


yAL 

0 

9/70 

13L/420 

0 13/35 



0 

-13L/420 

-L^/IAO 

0 -llL/210 

lVi05 

L 


..... (2.14) 

E IS the young's modulus of the material of the 


element, A is the uniform cross-sectional area of the link, 

L is length of the link and I is moment of inertia of the 
cross-sectional area, yAL mass of the element. The eflect 
of rotary inertia and shear are not included in the derivation 
of element matiices • The mass matrix given hero is a 
consistent mass matrix. 

2.3,2 SYSTEhi MODEL ; 

Section 2,2 describes an approach to represent 
the cam-f ollov/er pair into its equivalent lower pair system 
which can show the same dynamic and elastic behaviour as that 
of the actual system. Once this is accomplished the rest of 
the follov/er train along with this model can be represented 
by forming a system model to perform the defection analysis. 
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To construct the model of the entire system, system co-ordi- 
nates are to be defined at element node. (Fig 2.8) shows the 
system co-ordinates of the complete valve-train model. The 
system is consisted of total 7 elements and they are the 
numbers in the circles. It is noted tnat the system confi- 
guration that IS shovm in (Fig, 2.8) is for a particular 
geometric configuration of the cam profile which is taken 
care by the equivalent links. The system reference frame 
is 3(-Y, The crank inclines at an angle of 0 w.r.t, the 

The follower train has 5 elements. All are 
modeled as a planer frame element. The stiffness of return 
spring has its own spring index. Usually, the rocker arm 
has got peculiar shape. It is considered to be of two elements. 
One for the left rocker arm and the other for the right rocker 
arm. For modelling purposes they are approximated as beams 
with uniform rectangular cross-section. One can take more 
accurate model of the rocker arm by constructing a 2-dimensi- 
onal or 3- dimensional finite element mesh. After formulating 
a more precise model the no, of degree of freedom could be 
reduced [iS] so as to match element co-ordinates. Thus the 
basic analysis approach remains same but one pays heavy 
computational cost to obtain finer precision of the final 


answer 




INPUT CRANK MODELED 
AS CANTILEVER BEAM . 


IDEALIZED MODEL OF THE SYSTEM U,— -Ug 
SHOWS THE SYSTEM CO-ORDINATES. 


FIG. 2.8 
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The angle a { see Fig, 2,8 ) is a very small 
inclination of the rocker arm. In automative cams the 
total follower rigid body displacement is not high compared 
to the oimension of the links of the follower tram. The 
angle a actually takes care of the variation (small) in 
geometry. During rise or return period of the cam the rocker 
arm’s geometric configuration is specified by this angle a. 
Value of a can be obtained by dividing the rigid body dis- 
placement at tnat instant with the lower rocker arm length. 

Total of 9 system co-ordinates U is defined 
for the system model. To formulate the system equation of 
motions the element property matrices are to be generated 
in terms of system co-ordinates, so that the total system 
matrices are formed, 

. Fig, 2. 7 shows a generalised element in system 

oriented element co-ordinates. In case of a mechanism the 
different links take various geometric orientation throughout 
the course of the motion of input link. Thus element matrices 
are needed to be written down in svstem-coordinate form. 

The transformation from the element oriented element coordina- 
tes to system oriented element co-ordinates is accomplished 
by a transformation matrix [R], The general form of 


[R] Is [ ^ ° ] 

0 A 


,,, (2.15) 




ELEMENT 5 
(RIGHT ROCKER ARM) 


ELEMENT 6 (VALVE) 

FiG.2.9a 2.9f 

- SHOW NG THE INDIVIDUAL 

♦ ELEMENT CO-ORDINATES. 

2.9f 


Ci 
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where A ie given as 

SINCj) 0 

COS({) 0 where 4“ is as given in 

0 1 Fig. 2.7. 

The element mass and stiffness matrices for 
all the 7 elements are given in expressions (2.15a) to 
(2.15g) , These are oIotaDned for the corresponding 

element co-ordinates. The element co-ordinates are shown 
in( Fig 2.9a to 2.9f) for each element. Transforming them 
to system oriented element co-ordinates the transformation 
matrices are necessary. These transformation matrices [b] 
take into account the change in geometry. The six trans- 
formation jj^atrices ai® given in equations (2.15h “ 2,15m). 

The element coordinates and system coordinates 

are related by 

^u} = [B]^U[ ...(2.16,) 

for each element. 

The system oriented element mass and stiffness 
matrices are expressed as 

[k\ = [bJT [k]^[B]l ..,(2.17) 

[-J. = [B]J [m], [Bli ...( 2 . 18 ) 

i varies from 1 to no, of elements. 
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Element 


Element 


Element 


Element 


[m] = Ei.ll, to be calculated from (2.7) in 

section 2 . 2 . ^ (2.15a) 

[k], to be calculated from appendix II. 

1-1 

fk] ~ ^-1 1 here the coefficient are 

used considering longitudinal 
dffiflections. (2.15b) 

E^Ao 1 — 1 

[k] = -11 [_^ ^ ] (2.15c) 

3 

1/3 1/6 

[m] = YA 3 L 3 C 1/6 1/3 J 


[k] 


12E^l4 




6E4I4 




Symmetric 


«4A 


2^4 ^4 


4E4L4 


L 


'4 


...... (2.15d) 



Symmetric 


[m] = YA 4 L 4 


Element 5 


[m jrryA^L^ 420 


13/35 


Symmetric 


(2.156) 


Ck3 = 


^^ 5^5 


^%^5 


2E5I5 


I 2 E 5 L 5 Svmmetric 
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Element 6 


[m] 

[k] 


Element 7 


[m] 

[k] 


1/3 1/6 

Y ^ 6^6 ^ 1/6 1 / 3 ^ 


"6^6 



mass of return spring 
spring stiffness 


..,.(2a5f) 


(2,15g) 


since only one element co-ordinate is assigned to the spring. 


The transformations matrices are as follows: 


[b]j_ = {SinOji - 

CosO^ 

} 

.... (2.15h) 


CosOo 

SinOo 

0 

.... (2.15i) 

w 

11 

1 — 1 

z 


] 


0 

0 

Sin02 



1 

0 


.... (2.15j) 

[Bis = [ 

0 

] 

1 



CosO^ 

0 0 



[b]4 = 

; 0 

1 0 ■ 


.... (2.15k) 

1 

0 

0 1 




[B]5 

[Bis 


1 0 0 

0 CosO^O ■ 
0 0 1 


(2.151) 


(2.15m) 


0 


-1 
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It IS bo noted that [kj^and [m]j^ are 
function of input geometric parameter i.e. crank angle 0, 

Thus at every instant [k] and [m] are expected to vary. 

Having obtained the leleraent matrices in the 
global frame the next job is systematic assembly of all 
the element matrices so as to get total system mass and sti™ 
ffness matrices. 

The system mass matrix [M] and system stiffness 
matrix [K] are given by the following expressions 

NE 

[M] = 2 [5]. (2.19) 

1=1 ^ 

NE 

[K] = E [k]. ...,.(2.20) 

1=1 

V/here NE = r.j. of elements. The order of 
these matrices are equal to no, of system co-ordinate. 

Several computational schemes are available in order to 
assemble the element ma trices. e.g, , permutation vector 
method ^ code system of assembly [l3] etc. Here code-system 
is selected for assembly, because the flexibility in numbering 
the svstem cordinotes features it*s acvantage. 

Finally, v;ith this matrices determined the 
dynamic equation of motion for the system under stud/ is 
represented in the form as follows. 



^0 


= 

The equation (2,21) is written in terms unknovyn 
global CO— ordinates , It is actualy a set of differential 
eouations. The [C] matrix is system damping m^atrix, about 
which discussion is there in section 2.4,2, Vector ?f} 
contains the inertia forces due to the gross motion ( Rigid- 
body ) of the machanism plus any additional external forces 
if present. Equation (2,21) is the mathematical representation 
of the s’^stem model, 

2.4 A’lALYSIS APPROACH : 

2.4.1 METHODOLOGY ; 

The finate element displacement method is used 
efficiently in the olasto-civnamic analysis. A mechanism 
unlike a structure changes its position w.r.t, the fixed frame 
when it rotates. Thus the continuous rotation of the 
mechanism is divided in+o small time intervals. The continuous 
rotation can be more accurately taken care of in modelling 
if a co^Ginate is defined corresponding to the rigid bodv 

A 

degree of freedom of the mechanism. In the present work no 
such coordinate is defined at the rotation end and thus 
elealmination of rigid body degrees of freedom to make the 
system stiffness matrix non-singular is avoided. During 


l_Fj ( 2 . 21 ) 
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each interval the linkage lakes a certain geometric 
orientation according to its motion c laraoteristics . At 
each svch linkage configurata on the mechanism is converted 
into a structure by diminishing its rigid body degree of 
freedom. fee modelling of input crank as a cantilever is 
used to do this. 

At each such configuration the mass, stiffness 
orooerties of the mechanism treated as an elastic structural 
system are derived and assumed to remain unchanged during 
the interval. At this position the forces acting on this 
’structure’ are the rigid body inertia forces due to gross 
rigid bodv accelerations of the elements plus any external 
forces. The deflections are calculated solving the sot of 
O'-uations Oi motion. The next mechanism configuration is 
considered and the same procedure is followed. 

ElastO“dynamic analysis of a linkage using finite 
element approach essentially involves the following sequence 
of steps: 

(i) the idealization of the linkage structure is 
needed; this will require the selection of the type and 
size of the finite element to generate the system mesh. 
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(ii ) the system oriented mass and stiffness 
matrices are generated for each element. 

(ill) the element mass and stiffness matrices are 
superposed/ assembled systematically to develop mass and 
stiffness Oi the total structure of the svstem of linkage. 

(iv) determination of the unknown nodal displace- 
ments is made by solving a svstem of coupled ordinary differ- 
ential equations. 

(v ) all required numbers such as stress 
strains associated with the problem are computed. 

Previous section of this chapter (2.3.1, 
2.3.2) discuss about the implementation of first 3 steps 
given above. The solution procedure is discussed in the 
next section. 

2.4.2. SOLUTION PROCEDURE j 

The equation of motion developed in section 
2,3,2 is rewritten here, 

[M]fuj+ [C]|U}+ [K]|U} = ^Fj ...(2.22) 

This represent a set of coupled ordinary 
differential equations. The involving terms are already 
explained in previous sections. Damping effect is incorporated 
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in tlie j orn of s'^ctem damping matxTx [C], The complete solution 
oi the n'uai.aon i'1,2'^) consists of accomplishing 3 parts. 

(i) determination of eigenvalues and eigenvectors 
o' tho coi responding homogeneous equation. 

[i.0m+ 0 (2.23) 

T;v-> of fuel o: damping on the eigenvalues and eigenvectors 
is nssu:"^d to bvi negligible. 

(ii) transformation of this set of coupled equ- 
ations into a sot of uncoupled equations expressed in the 
form of normal co-orc'inates. 

(iii) solve for the complete solution of the vector U. 

Tho implementation of the first step that is 
tho cic' lormin-.t ion of eigonsvstem is described in Appendix III. 

To trsnsform the set of couoled equations into 
a set of uocouolec’ eouations the svs tern co-ordinates U are 
convort'^c’ to nornnl coordinates as shown below 

fu] = [(t)]^nS ....(2.24) 

whoro {r) { t )} 1 s the displacement vector in normal coordinate. 

[cj)] represents the modal matrix which contains resulting 
eigenvectors columnwise as calculated in step (i)* [({)] 

remain constant during the choosen interval. 
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Successive differentiation of equabion (2.24) 

yiolcis i/Hc r oilowin.Q’ "two OQUB'txons 

u = [<|)] ....(2.25) 

and U = [(|i] |7||. ....(2.26) 

suos ti tu ting equations (2.24), (2,25), and (2,26) in equation 
(2,22) it can be rewritten as 


[M] [(tl]Ylij+ [c] [K] [(t)]Jn} =-fFj ....(2.27) 

Premui Liplying this equation by and using the orthogona- 

litiGs properties the equation formed is as in (2.28), 


where, 



[M] = [(!>]’■ [M] Wl 

[K] = [K] [(])] 

[5] = [<t>f [C] [((>] 


(2.28) 


... (2.29) 


All the matrices [M], [K], [C] are diagonal matrices. Generally, 
the system damping matrix [C] is taken as proportional to 
either [m] or [k] or a linear combination of them. Obviously, 
the matrix [C] will also to diagonal. But the use of a single 
constant to represent the damping characteristics of the 
systems with multiple degrees of freedom is not realistic. 

For this reason, the standard practice is to consider [C] as 
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oroDor Lionr"’! 1o tho crita cal damping of the system at the 
nor»'’nl jnoi'G [9]. Again the matrix [C] becomes diagonal and 
n'''’JtTon (2,29) becomes uncoupled. These uncoupled equation 
assume thf^ '"ollowing form, 

m.Ti^(t) + 2m.^. m^Ti^(t) + = p^...(2*30) 

i = l,2,...m, ; ffi = no, of modes foundout in 

oioen system. 


The above equation represents the i-th equation 
of tho set of resulting uncoupled differential equations. 

The associated terras are 


% 

^1 


P 

i 


Cl). 

1 


2m. CD. 
1 1 



= displacement at the i-th normal coordinate 
= velocity at the i-th normal coordinate 
= accelerations at the i-th normal coordinate 
= i"th diagonal of [M] 

= i-th diagonal of [k] 

= i-th element of 


i-th natural frequency 

critical damping at the i-th normal 

damping ratio (critlcaT'dtSpilg' 
normal mode , 


mode 
the ith 
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i.iatrix [Cj for the choice of daraping in (2.30) is 


[C] = 2 




^ 2^2 


w V 
m > m 


[m] ...(2.31) 


As the damping is choosen to be proportional to the critical 
damping the iAai,rices [c] and fc] are never required in 
actual computation. 


Noting that [K] =[w^] R] ...(2,32) 

whore [w J is also a diagonal matrix , equation (2,30) can 
be rewTitten as 

’nh(t) + 2a)^'^^n-^_(t) + il^(t) = p^ ...(2,33) 

^ P. 

where, p^ = ...(2,34) 

Each element of p^ is assumed to remain contant during the 
current interval. 


The solution of the equation ( 2.33 ) is given as [ 17 ] 


ri3,(t) 




e 




^ [ A, Cos(ii. t) + 


y— Sin (^j^t)] 


7i 


....(2.35) 

Where first term is the particular Integral, and 
second term is the homogeneous solution. 
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. . . . (2.36) 
.... (2.37) 

.... (2.38) 


11^(0) and rj^(O) are the ith initial value of the 
i sp] ■•cenenL end velocity of the vector t). 


The displaceiient T] determined from equation ( 2.35) 
are crniioformed back to svstem coordinate U with the aid of 
tlK' oinjetion (2.24). Finally the unknown nodal displacements 
r'iTO cnlci'lated out. If one wishes he can also find out 
the r.v£ Lem oriented element co-ordinates employing rhe 
matrix used for assembJv of all element matrics. 


2.'^. 3 rOiiCE vector ; 

The force vector CF^in enuation (^.22) consists of 
inertia forces due to rigid body motion. {F^ is obtained 
by multiplying the system mass matrix with negative of 
acceleration vector . 

= - [m] ... (2,39) 

where , R = rigid body acceleration vector along the 
system coordinate. In this approach variation of acceleration 




over ti.c link lengths is not considered. {fI during each 
iniorv-‘J is ''ssurned to renain constant. 
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It can be said from section 2,1 that the 
acc-^J erati on of the roller centre ( see Fig. 2.1 ) and 

that of the slider are same instantaneously. The slider 
acceleration can be computed from the expression ( 2 , 40 )USj- 

o _ ^ Cos ( J 80—202 ) +Cos 

3 rr — . (j) 2!^ , Sin 0Q “h 

(n2-Cos%)=/2 

. ... (2.40) 

whore, = constant angular velocity of the crank 



a = slider ac ileration . 
s 

Acceleration of the push rod is same as the 
roller centre. The valve-end of rocker arm moves with an 
acceleration. 

a=-CXa (2.41) 

V s 

where, C is the lever ratio of rocker arm,’-* sign indicates 
the opoosite direction. 
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/t this luction the return spring force is also 
.■‘ct'ive ■'hie’'' Ts the pror’ucl of soring stiffness and spring 
r’ 'flee, ion ( rigid ). This force is added with corresponding 
olorro'''! in the vector |f} and the resulting {fJ- is obtained. 
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CHAPTER HI 
CO-'IPUTCR PROGRAi'-iS 


For the purpose of implementing the steps 
pr''‘S<'r.t in tne section analysis approach a 

cohinutor program has been c’evelopec’. The aim is to find 
out 1 m s >li't3on of the euqations (^,22 ) which in turn 
produces the displacenent vector The program is 

i"riltci in fORTRAM-IV and run on DEC~1090 system. 

BASIC CC.TUTHIONAL STEPS; 

(i) Read all the data describing the geo- 
metry Oi the follower linkage ( i.e. link dimensions), 
matorinl property , rotur' spring stiffness etc. The data 
related to element assembly are also read in input. 

(ii) Get the data for individual masses and 
lengths of the equivalent links of the cam-follower through 
a seperate subroutine. 

(iii) Compute the stiffness and mass properties 
of each element employing equations (2.15 a - 2.15 g ). 

(iv) Starting from an initial position of 
the cam ( the initial position or the crank is derived from 
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’• *' ! 1' orn une Glorient matrics in step (iii) 

tc - ‘I'" ' '-''''.toe’ co-orc'inato system. 

(v) Asaenblo all the system oriented element 
'' ’ '■'i lotal system mass and stifness matrices 

[ ',j irr. CKj respectively. 

(vi) Solve the eigen value problem as described 
in p viion ('^,23) to find all the e3gen“values and find the 
",Ddnl .1 trix [cj)] consisting of eigen vectors, 

(vii) Compute the force vector (section 2,4.3) 
in global co-ordinate frame, 

(viii) Uncouple the svstem governing equations 
(2.22) and solve for the displacement vector g in normal 
co-ordi .*‘toG. 

(lx) Transform the vector to obtain by 

the expression (2.24). 

(x) Go to the next cam position by adding the 
incrernonLal angle with the current angle and repeat the 
entire stops from step (ii) to step (x) except siep (iii) 
unta 1 the final cam position is reached. 
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SUBROUTINES : 


Followinq subroutines are used to back the 

main program. 


(1) ELivICON : This subroutine forms the element 

mass and stiffness matrices in element 
oriented co-ordinate system and then 
convert them to system oriented frame. 
This routine is invoked by the routine 
SYSiJiff. 

(2) TRANS : This routine is used to calculate the 

transformation matrix corresponding to 
each element for successive cam position 
or crrik position. The routine ELMCON 
calls this routine, for transformation, 

(3) SYS' lAT ; This subroutine servos the purpose of 

systematic assembly of all element 
matrices provided by routine ELMCON. 


(4) JACOBI; This routine solves the generalised 

eigenvalue problem, AX=)‘BX . It makes 
use of the generalised Jacobin algorithm 
The output of tha s routine is modal 
matrix [(j)] of the order NSCXNSC where 
NSC IS total no. of svstem co~ordinate 
and eigen vector > of order NSC, 
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(5) SOLVn : "ii'' Golulion of the System^s dynamic 

cojjutjon of motions (2*22) is done by 
tiiis ToutDno, It invokes the routine 
JACOBI and chock for the orthogonalities 
properties i,G,^ wf m [(f] = [>■] 

and [([i]'’' [m] [(j)] = [l], [l] Is an 
identity matrix , The vectors *0 and 
final displacement vector U are then 
determined by this routine. 

(6) EQLINl; The equivalent link lengths and their 

included angle v;ith reference to 
fixed frame are produced by this 
routine. The rigid body displacement 
a^ thf^ cam end is also output of 
this xO\>tine, 

(7) FOjISV: fh'.s subroutine is used to calculate 

force vector representing the R,H,S 
of th''' equation of motion (2»22), 

Apart fro'i the above routines function 
subprograms are used to determine stiifness of the crank, 
connecting rod in the equivalent mechanism and mass of the 


crank* 
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CHAPTER IV 
COiVPUTATIOMAL WORK 

To judge the workability of tho analvtical 
treatnents presented in Chapter II a numerical example is 
solved by means of a computer program ( Cnapfcer III). 
Computational observation and discussions are depicted in 
this chapter. 


A.l NUraUCAL EXAMPLE; 


Tho numerical example that is taken up contains 
a cam of the tvpe DRRD ( i.e., Dwell-Rise-Re turn-Dwell ) which 
IS verv common in automotive overhead valve train linkage. 

Tho can is shown in (Fig. ".1 ). It is comprised of 3 circular 
arcs for the rise-return portion. Other input data pertaining 
to tiio rest of the mechanism that is of the follower train 
are given below : 


Length of tho Pushrod - 22.0 cm. 

Diameter of the Pushrod — 1,5 cm. 

Length oF the left roker arm - 4,5 cm. 

Cross Sectional area - (2.0 x 1.0 ) cm' 

Length of the right rocker arm- 5.5 cm . 

Cross -section area “ (2.0 x l.O) cm , 

Length of the valve “ 9.0 cm 

Dia. of tho valve “ 1*5 cm 
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Return spring stiffness 

Roller radius 

Young’s itiodulus of Steel 

Poissons ratio 

Material density 

Damping ratio is taken es 


~ 40.0 kgf/cm, 

- 0,8 cm, 

~ 2,1 X 10^ kgf/cm^, 

- 0.3 

~ 7,8 X lO”^ kg/cm^ 

«- 0.05 


The required data’s for the earn profile and follower are 
shown in Fig. ^’.1, 


4.2 RESULTS AND DISCUSSIOI^Sj 

The cam profile that is being considered 

in this example shows 180 degrees of dwell characteristics. 

Durina the dwell period the follower does not execute any 

motion. Thi's the analysis is aerformed w thin a range of 

180 degree. The start of the rise ( point in Fig. 4,1) 

initial 

IS considered as^osition of cam ( 0 degree ) and the end of 
return ( point ) as final cam position (180 degree ), 

centre of curvature of respective 
profile region. But actually in one full rotation of cam 
the initial position is taken from dwell followed by rise, 
return and again dwell. For sake of analyses the cam rota- 
tion variation is illustrated from 0-180 degree in all 
figures. The analvsis is carried out with an interval of 
5 degree. The linkage js initialised for zero displacement 
velocity. 




CAM ROTATION (DEGREES 
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Tho moss and sti'Ffness values of the crank 
in the equivalent linkage for the present example are 
computed using the formula (2*7) and (2,9). ’''hose values 
are 0,5428 Kg and in the order of 10^ kgf/cm respectively* 
The stifiness assigned to the connecting rod is found (eq, 
2,12) out to be 0,3666 x 10^ kg/cm* Though the connecting 
rod ( link2 )is assumed massless to avoid computational 
instability very small value of mass is consic'ored for it. 

The results of the deflection anal'''’sis are 
plotted for the node 9 (Fig. 2,8) that is at the valve end, 
Oftonlv for automotive valve train cam linkage the displace- 
ment of the end mass 3 s of paramount concern. Adherence to 
a prescribed motion characteristics of the end element i.e,, 
of the valve determines tae proper functioning of the entire 
unit. The error ( or the difference between the cam command 
and follower response ) may be represented by deflection 
at node 9 i.e., Ug. Added with the rigid body displacement 
the net value may represent the total displacement of valve 
end . 

The displacement characteristic of the cam 
profile is shown in (Fig. 4,2 ). The figure represents the 
rigid body displacement at the valve end for the rise-return 
portion. Fig. 4,3 shoves the characteristics curves of 
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the circular are cam one translating roller folloVi/er for 
the example being considered, Tho displacement and accele- 
ration are plotted against the cam rotation angle at tho 
cam end. The curves are shov'n for both tho rise return port- 
ion and the rise is started to be shown from 0 degree. 

But actually in the case of a DRRD cam this would be 90 
and the end of return 270*^, It is evident from tho figure 
th^t the accleration suddenlv changes from positive value 
to negative value and the maximum negative acceleration 
is greater than the maximum positive one. This is obser- 
v>sd at the blending point Kj^ and at (K 2 ) » because at this 
particular point two equivalent slider-crank mechanism 
represent tho cam motion. The symmetry is notable because 
of the svmmetric shape of cam profile. 


The elastic deflection at valve end (U^) obt- 
ained by solving equation (2,21) and computed for each 
interval is shown in(Fig ^.4 and 4.5), This i$ plotted 
against the cam rotation and for various operating speed 
of the cam. At the starting of rise ( 0 deg. ) the dof lectio 
■j c as zero. Few observation that are moce are as 


follows : 


(i) The shape and nature of the deflection curves 
resembles verv much to that of acceleration curve (Fig. 4.2) 
This is expected because the forces ( inertia force), that 
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cause deflections vary in the same manner as ‘chat of the 
acceleration curve against various successive geometric 
configuration, 

(li) The elastic deflection \Jg appears to remain 
almost constant on the portion and K 2 K 2 for various 

input speed. This deflection assumes higher value but in 
the opposite direction in nose portion (K 2 ^K 2 ) of the cam 
compared to the other region of the profile. The marked 
differences of the deflection value on the two region may 
be attributed to the nature of acceleration and consequently 
the force in this region of cam profile . Equivalent link 
lengths are different that is the instantaneous ’structure’ 
is also changed which changes the system stiffness and 
mass matrix. The difference between the maximum deflection 
and minimum deflection in the region l^i^2 ^ 
more pronounced at high speed (1800 rpm) compared to that 
at low speed (800 rpm). The curve shape is more flattened 
as the speed decreases. From the obtained result (not 

given ) it is seen that deflection at low rpm (500) becomes 
too small and thetotal displacement curve at valve end almost 
coincides vath the kinematic displacement curve. 

Fig. 4.6 describes the deflection at node 9. 
These deflections are based on static analysis with the 
rigid body inertia forces as the force vector.- No Eigenvalue 




- 0.8 
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problem is solved. Those deflection are seen to be smaller 
than that solved through dynamic equation. 

The natural frequencies of the linkage depend 
on its instantaneous configuration. The varying geometric 
configuration causes the natural frequency also to change. 
The nature of variation of first natural frequency with 
different cam position is illustrated in the (Fig. d.?). The 
variation is seen to remain almost constant, in the flank 
portion, while variation is more prominent in nose region. 
This indicates that the change in configuration of the 
input crank afreets insignificantly the system mass and 
stiffness prooerties in the region and K2K2. 
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CHAPTER V 
CONCLUSIONS 

5.1 TECHNICAL sUMlvlARY : 

The tochnique of kinetoelastodynamic 
analysis of niechanism finds its importance in a situation 

where a mechani om ,• x , ^ . 

^ism IS operated at high speed. The accuracy 

requirem nts ij-j output displacements, considerations of 

elasticity and dynamic factors in producing stress, strains 

deflection in high speed planar mechanisms are major problem 

which the rigig^^^^ assumption cannot overcome, KED, an 

advanced analysis and synthesis procedure is efficient 

in handling a circumstances which takes into accounts links 

flexibility, high inertia forces and provide means to 

obtain accurate solutions. 

The applicability of the method was mainly 
limited to lower-pair planar mechanisms. Higher pair 
mechanism remain neglected in the context of using the 
method for better analysis. Cam mecnanism, fall under 
higher pair classification and stand as an important part 
in the construction of modern day machinery. Evidences are 
plenty regarding the use of the mechanism. They are seen 
in automotive engines, aircraft engines, textile machines 
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and so on, A cam~actuated valve train linkage found in 
every automotive engine plays a responsible role in assisting 
to achieve the smooth operation of the engine. This, on the 
other hand, largely depends on the correct synthesis of the 
cam profile. The discrepancy existing between cam command 
and follower masses owing to the flexibility in the linkage 
needs correct evaluation for the betterment of the cam 
performance. 


Conventional methods using lumped masses, 
spring etc, fail to provide a better analysis because of 
their certain inherent drawbacks. The necessity arises to 
exploit the potentiality of the newer methods so that a 
satisfactory and good analysis may be performed for cam acti- 
vated valve train linkage. One faces difficulties while 
formulating an elastodynamic model of a cam mechanism because 
of the fact that the cam along with the follovi/er does not 
match with the typical element shapes in elastodynamic liter- 
ature in forming the cam as a known element. Certain 
considerations must be accounted so as to make the correspon- 
ding model of the cam follower possible and one can go ahead 
to model the entire linkage and perform the elastodynamic 
analysis. Thus some attempts may be tried in carrying out 
the analysis for the estimation of elastic deflections. 
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The obiective in the present work is achieved 
in tne following 2 phases ; 

(l) The equivalent piodel of the cam follower pair 
is proposed Vi'hich shows kinematic ecuivalencv and produces 
the approximatelv same elastic and dynamic behaviour as 
that of the actual svstepi. The stiffness and masses of the 
links are derived from the above considerations and values 
are imparted to the equivalent links. These masses and 
stiffness do contribute towards the formation of total 
SYS tern mass and stiffness matrix. 

(2) The elastic model of the entire svstem in 
terns o" global co-ordinatos is created. The equation of 
motion of the svstera is solved for the unkown nodal displace- 
ments. These displacements represent the elastic deflection 
at corresponding nodes. 

The findings of the present work may be summa- 
rised as follows; 

The elastic deflection at the valve end (error) 
increases with the increase in input speed. The nature of 
deflections depend on the cam profile characteristics. The 
natural frequency of the instantaneous structure does not 
vary largely with the variation of inout crank position. 
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RECOMj-ieijdatIONS; 


In the beqining chapter the eraphasis is placed 
he fact that the present v'ork is an attempt to analyse 
- lOo toe ynamically a cam-modulated linkage. In view of the 
Ifiiitationslisted in section 1,5 , the input crank is modeled 
fis a cantilever beam. The modifications can be made by 
incorporating a coordinate at the suoport end of the crank 
thus taking care of the rigid bodv d,o,f, and the mechanism 
can be analvsed under actual boundary conditions. 

The ocui violent model to replace the cam™ 
follower pair can also oe refined for more accurate repres- 
entations of the dynamic and elastic behaviour of the actual 
cam-follower. Consequently the raass and stiffness assigned 
to crank and connecting rod will be more improved. 

No distribution of the rigid bodv accelerations 
over the element length is assumed. Instead the inertia 
force ( or body force ) can be derived at the element level 
for each element and cam be transformed to system reference 
frame to obtain the final force vector by assembling. 

In this vi/ork each link is regarded as an 
element. For more accurate modelling particularly the pushrod 
and the valve can be subdivided into few elements inspite 
of one. 
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Appendix I 

Determination of area moment of inelPtla of cam profile 


A cam profile is shown in Fl0» (A*- 1 ) which is 
pivoted at the point 0. This profile is a convesi profile* 
As mentioned in equa tion(?, 3 ) the polar moment of inertia 
of the area ( Jc,area ) is to be computed* 


Consider a small sector bounded bv r^ and 
The included angle is aq:* The polar moment of inertia of 
this small sector mav be given as Aa “4 about 0 


where r^ s: and sector is 

assumed as a small circular sector with mean radius r^^. The 
total moment of inertia V'ill be given [s] bv 


To calculate I the values r and 9 are required 

for n no. of points on the profile. When a cam profile is 

known completely r 9. can be computed for all i. For a 

i ^ 

cam with dwell, only in the rise-return portion the values 
of r,0 have to be calculated. 
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Appendix II 

Derivation of Cam Stiffness 

The derivation of approximate estimate of 
stiffness of the cam bodv is given here. The(Fig 2, 3^2. 4) 
shovjs a earn profile and a tapered beam which approximates 
the cam profile. This beam has its cross sectional area 
varying over the length of the beam. Let’s say a load P 
is acting at the tip. Under this load the maximum deflection 
can be computed using castigliano ’s theorem [n]. 


d = 

where d= 
U= 
L= 
t= 


L p p 

dU r P X dx ... 

deflection under the load P 
strain energy due to P 
length of the beam 
uniform thickness of the cam (beam 


) 


(A21) 


X is measured from the tip towards the support 


L 


=P 


/ 


2 ^ 

X dx 

E. ^ ( a+cx) • 8 


.. . . (A22) 


where C= 


L 
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or 


d = 


I .E.t/^ [ 1„(|) - (1- I-) + |( 1- % )1 

( A23) 

Thus the stiffness from(A23) can be 


computed to be 


S' E t 

K. 2 J ,..<A24) 

[ in (|) - (1- |) . I (1- ^)] 

This is an approximate estimate of stiffness 
actual cam which will give a feel of the order of magnitude 
of the value of Ky,. Kj^ is dependent on the dimension, b, a 
and L. These values are to be selected compatibly with the 
size of cam , 



75 


Appendix III 

Determination of Eigen-value artd Eigen'-vector* 

The undamped equation of motion given by 
equation (2.23) is a generalised eigen-value problem, in 
the form, 

K (j) ^ >M ({) .,...(A31) 

where , 

K is a (nxn) stiffness matrix 

M is a (nx) mass matrix 

A 

(]) is a (mxn) modal matrix containing 

m no, of eigen vectors corresponding to 
m eigon values, 

is (mxm) diagonal matrix whose non-zero 

terms are square of natural frequencies. 

matrix 

In the present problem the order of^(i,e, n ) in 
K and M are comparativelv small. Hence the modal matrix 
(j) may contain all the eigen vectors and equation (A31) can 
be solved for all eigen values. 

There are many methods for solving generalised 
eigen-value problem ( e,g, generalised Jacobi method, Q-Z 
algorithm etc.). The generalised Jacobi method is suitable 
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v'nen the order of the problem is small. Thus this method 
is choosen for solving eq (A31), Besides, it has other 
advantages [l] like it can handle ill-conditioned matrices, 
it avoids solution of standard eigen value problem and 
solves (A31) directly, fast convergence is expected when 
the off-diagonal elements are small. 

The modal matrix (}) satisfies the M-'^rtho- 
gon^^lities as follows ({)^ MC]) = I - (A32) ,.,.(A32) 

and ({)’^K(j)=’X .,..(A33) 

where I is an identity matrix (A32) and A33) can be used as 
a cheek of the solution. 
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